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Abstract. We present CCD photometry and a proper motion 
study of the young open star cluster, NGC 581 (M 103). Fitting 
isochrones to the colour magnitude diagram, we found an age 
of 16 ± 4 Myr and a distance of roughly 3 kpc for this cluster. 
The proper motion study identifies 77 stars of V — 14.5 mag 
or brighter to be cluster members. We combine membership 
determination by proper motions and statistical field star sub- 
traction to derive the IMF of the cluster and find a quite steep 
slope of r = -1.80. 

Key words: open clusters and associations: individual: 
NGC 581 (M103) - astrometry - stars: kinematics - 
Hertzsprung-Russell and C-M diagrams - stars: luminosity 
function, mass function 
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1. Introduction 

The shape of the initial mass function (IMF) is an important 
parameter to understand the fragmentation of molecular clouds 
and therefore the formation and development of stellar systems. 
Besides studies of the Solar neighbourhood (Salpeter 1955, 
Tsujimoto et al. |1997 ), work on star clusters plays an impor- 
tant role (Scalo 1986| ), as age, metalicity, and distance of all 
cluster stars can be assumed to be equal. 

Most of the previous studies indicate that the IMF of a star 
cluster has the shape of power laws 

N{m) - (1) 

within different mass intervals. Th e foll owing typical values of 
their exponents as given in Scalo (1998): 



r = -1.3 for TO > IOM0, 

r = -1.7 for IMq < to < IOA/q, and 

r = -0.2 for TO < IMq. 



(2) 



Knowledge of membership is essential to derive the IMF espe- 
cially of open star clusters, where the contamination of the data 
with field stars is a major problem. Two methods for member- 
ship determination are in use nowadays and each of them has 
its advantages and disadvantages: 

Send ojfprint requests to: Jorg Sanner, jsanner@astro.uni-bonn.de 




Fig. 1. 60 s y filter exposure of NGC 581 (M 103) taken with 
the Im Cassegrain telescope of Hoher List observatory. The 
field of view reproduced here is approx. 15' x 15'. North is 
up and east to the left. The stars marked as A, B, and C are 
mentioned in Sect. 4.2 



- The classical method is to separate between cluster and 
field stars by their proper motions: All cluster stars can be 
expected to move in the same way, whereas the field stars 
show more widely spread and differently centred proper 



motions (see e.g. therecent workof Francic 1989). For each 



star a membership probability can be specified. To obtain 
a sufficient epoch difference, old as well as recent photo- 
graphic plates are needed to measure proper motions, so 
that this method is limited by the comparably poor sensi- 
tivity of the old plates. 

With the introduction of CCD imaging to astronomy, statis- 
tical field star subtraction became more popular Assuming 
(almost) identical field star distributions in the cluster re- 
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gion itself and the surrounding area, the distribution of the 
field stars can be subtracted from the one of the (contam- 
inated) cluster area. This makes sense for fainter stars but 
with the bright stars one deals with statistics of small num- 
bers. 

Our work combines these two methods of membership de- 
termination: The proper motions are investigated for the bright 
stars of the cluster, whereas the fainter stars are treated with sta- 
tistical field star subtraction. From the cleaned data we derive 
the luminosity and mass functions of the cluster. 

NGC 581 (M 103), which is located at asooo.o = l''33.2™, 
^2000.0 = +60°42', was chosen as a first test object for our 
technique because the cluster and a sufficiently large field star 
region can be covered within the field of view of the telescope 
used. A V image of NGC 581 is shown in Fig. |l]. 

In Sect. H we present our CCD photometry, and in Sect. ^ 
a proper motion study of NGC 581 and another cluster which 
is located on the photographic plates, Trumpler 1 . The result- 
ing colour magnitude diagram (CMD) is discussed in Sect. ^ 
leading to the determination of the IMF of NGC 581 in Sect. 
141 



2. CCD photometry 

The photometry is based on 22 CCD frames taken in Johnson 
B and V filters at the Im Cassegrain telescope of Hoher List 
Observatory. The telescope was equipped with a focal reducing 
system and a 2k x 2k CCD camera called HoLiCam (Sanner 
et al. 1998 1, which has a pixel size of 15/xm x 15/im and a 
resolution of 0.8"pix^^. The field of view covered in this con- 
figuration is a circular area with a diameter of 28'. Information 
about the images used for the photometry are summed up in 
Tab. [l|. 

The images of equal exposure times were averaged, result- 
ing in integrated exposure times of 35 min in V and 60 min in 
B for the longest exposures. The shorter images were used to 
gain information about the bright stars which were saturated af- 
ter longer exposure times. After standard image processi ng the 
photometry was performed with DAOPHOT II (Stetson 1991) 
running under IRAF. After an error selection process, the data 
were calibrated from instrumental to John son m agnitudes using 
the photoelectric sequence of Hoag et al. ( 1961 ). Their standard 
stars as well as our instrumental magnitude values and their de- 
viations are given in Table |[ 

We applied the following equations to transform instrumen- 
tal to apparent magnitudes: 



v-V ^ aQ-ai-{B -V) 



{b-v)^{B- V) 
with 



a'l -{B-V) 



ao = 5.536 ±0.04, oi = 0.094 ± 0.02 



= 2.418 ±0.04, 



= 0.139 ±0.02 



(3) 
(4) 

(5) 
(6) 



where B and V represent apparent and b and v instrumental 
magnitudes, respectively. Mean photometric errors in different 



Table 1. Summary of the CCD images from the Im Cassegrain 
telescope at Hoher List Observatory used for the photometry 
and the proper motion study 



filter t 



exp 

[s] 



number of exposures 
photometry proper motions 



V 
V 
V 
B 
B 
B 



10 
60 

300 
10 
60 

600 



Table 2. Internal standard stars of Hoag et al. (1961 1 with the 
deviations of the computed from the catalogue magnitudes. 
Stars 1 and 2 were saturated even in the shortest exposures 



no. 


V 


B-V 


AV 


A(B- V) 




[mag] 


[mag] 


[mag] 


[mag] 


3 


9.09 


+0.21 


-0 


020 


+0.039 


4 


10.45 


+0.24 


-0 


038 


+0.007 


5 


10.59 


+0.17 


-0 


013 


+0.002 


6 


10.81 


+ 1.93 


+0 


004 


-0.034 


7 


11.22 


+0.23 


-0 


028 


+0.015 


8 


11.35 


+0.26 


-0 


002 


+0.013 


9 


11.76 


+0.22 


-0 


016 


-0.022 


10 


11.84 


+0.25 


-0 


031 


+0.021 


11 


12.34 


+0.10 


+0 


083 


-0.105 


12 


12.76 


+0.24 


+0 


090 


-0.047 


13 


13.17 


+0.30 


-0 


016 


+0.023 


14 


13.21 


+0.26 


+0 


012 


-0.119 


15 


13.27 


+0.44 


+0 


034 


+0.019 


16 


13.46 


+0.30 


-0 


009 


+0.031 


17 


13.47 


+0.35 


-0 


028 


+0.075 


18 


13.59 


+0.77 


-0 


001 


+0.060 


19 


13.76 


+0.30 


-0 


008 


-0.010 


standard deviations 


0.036 


0.051 



Table 3. Photometric errors in different magnitude ranges 



range 



Al/[mag] AB[mag] 



V< 12 mag 0.007 0.009 
12 mag <V< 16 mag 0.012 0.021 
16mag<l^ 0.028 0.092 



magnitude intervals are given in Table ^. One can see that the 
errors for B increase more rapidly as a consequence of HoLi- 
Cam's poorer sensitivity in blue wavelengths. Although the to- 
tal exposure time in B was almost twice as large as in V, the 
limiting magnitude of the photometry is defined by the B im- 
ages. 

From these data, we determined a CMD which is shown in 
Fig. ^. It represents a total of 2134 stars for which both V and 
B — V magnitudes are available. We present the photometric 
data of all these objects in Table Q 

The CMD shows two main sequence features and a scat- 
tered giant branch in a colour range around B — V ^ 1.4 mag. 
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Table 4. List of the photometric data of all stars measured in 
the CCD field of NGC581. Star numbers higher than 10,000 
account for objects for which proper motions were determined, 
too. Only the ten brightest stars for which both photometry and 
proper motions are available are listed here, the complete table 
is available online at the CDS Strasbourg archive 



Table 5. Photographic plates from the Bonn Doppelrefraktor 
used for the proper motion study. The plates marked with MS 
were digitized with the PDS machines in Miinster, the ones 
with T with TPS in Tautenburg. For comparison, both digiti- 
zations of R 0285 (Miinster and Tautenburg) were used during 
data reduction 



6 
7 

8 

9 
10 

11 

DO 

13 
14 
15 
16 
17 
18 



No. 


X 


y 


V 


B-V 


plate no. 


date 


iexp 


scan 








[mag] 


[mag] 






[mill] 




14380 


1024.178 


1122.116 


10.603 


0.168 


R0285 


23.01.1917 


30 


MS/T 


14561 


1136.343 


1327.688 


10.806 


1.964 


R0294 


09.02.1917 


30 


MS 


14503 


1099.685 


1056.475 


11.028 


0.177 


R0295 


09.02.1917 


3 


T 


14411 


1046.961 


1088.074 


11.432 


0.179 


R 1281 


03.12.1977 


60 


MS 


14288 


972.818 


1087.454 


11.478 


0.191 


R 1291 


05.12.1977 


64 


MS 


13899 


697.822 


446.703 


11.674 


0.530 


R 1295 


19.12.1977 


4 


T 


13894 


708.770 


1589.465 


11.776 


1.917 


R 1298 


19.12.1977 


4 


T 


14223 


925.276 


1073.279 


11.776 


0.242 


R 1899 


24.03.1998 


60 


T 


14341 


1003.143 


1165.883 


11.785 


0.170 










14928 


1391.922 


1103.418 


11.871 


0.229 












1 1.5 
B-V [mag] 



Fig. 2. Colour magnitude diagram of all objects detected on the 
CCD frames of NGC581. Note the two main sequence struc- 
tures and the scattered red giant branch. The origin of these 
features is discussed in Sect. Q 

More detailed analysis of the CMD is being presented in Sect. 



/ = 5.1 m instrument represent a region of 1.6° x 1.6°. The 
plates were digitized with the PDS machines of the Astronomi- 
sches Institut Miinster and the Tautenburg Plate Scanner (TPS) 
of Thiiring er La ndessternwarte Tautenburg (Brunzendorf & 
Meusinger |1998| ). The positions gained with DAOPHOT II 
from 15 CCD frames were added to the plate data. Tables [l| 
and H give an overview of the material included in the proper 
motion study. 

The celestial positions of the stars were determined from 
the plate coordinates with respect to six HIPPARCOS stars 



(ESA 1997 1 using an astrometric software package devel- 



oped by Geffert et al. (1997). We obtained good results us- 
ing quadratic polynomials (6 plate constants) for transforming 
(x, y) to {a, S) for the photographic plates and cubic polynomi- 
als (10 constants) for the CCD images, respectively. The mean 
positional deviations of the HIPPARCOS stars after the first re- 
duction step were of the order of 0.1" in both right ascension 
and declination. Using the output positions and proper motions 
of each step as the basis of the next run, we derived a stable so- 
lution of proper motions of a total of 2,387 stars on the whole 
field after four iterations with a mean error for the proper mo- 
tions of approx. 1.1 masyr^^ in both coordinates. The differ- 
ences between our measurements and the HIPPARCOS data 
("observed - calculated" or O — C values) are listed in Table |[ 
Compared with our measurements, HIPPARCOS star no. 7155 
showed high deviations — probably caused by a double star 
nature of this object (Wielen et al. 1998) — and was therefore 
excluded before the data reduction. We present the proper mo- 
tions of all stars in Table ^. 



3. Proper motion study 

3.1. Data reduction 

For the proper motion study eight photographic plates from the 
Bonn Doppelrefraktor (until 1965 located in Bonn, thereafter at 
Hoher List Observatory) were used, covering an epoch differ- 
ence of 81 years. The 16 cm x 16 cm plates of the D — 0.3 m. 



3.2. NGC581 

228 stars were located within 10' from the centre of NGC 581. 
Only these shown in Fig. ^ were taken into account for the vec- 
tor point plot diagram. Due to the cluster's large distance (see 



Sect. X.2), the centres of the distributions of field and cluster 
stars are more or less the same. Therefore, the separation be- 
tween field and cluster stars is not very apparent. 
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Table 6. "O — C" values of the measurements of the HIPPAR- 
COS stars which were used for the transformation from pixel 
to celestial coordinates 



star no. 


Aa 


AS 


A^a COS 5 


Afj.s 




["1 


["] 


[mas yr~^] 


[mas yr~^] 


6979 


0.027 


-0.020 


0.22 


-0.95 


7640 


-0.006 


-0.080 


-2.20 


-1.38 


7497 


-0.006 


-0.130 


-0.91 


-1.72 


6793 


-0.003 


-0.020 


1.17 


2.40 


6927 


-0.011 


-0.080 


1.93 


-1.10 


6841 


0.050 


-0.330 


-0.85 


0.50 



Table 7. List of all proper motions determined from the pho- 
tographic plates of NGC 581. The positions are given for the 
epoch 1950.0 in the equinox 2000.0 coordinate stystem. The 
stellar numbers are the same as in Table 0. Only the proper mo- 
tions of the same stars as in Table Q are presented here, the com- 
plete table is available online at the CDS Strasbourg archive 



No. 


a 


S 


fla COS 3 


f^6 




[h m s] 


[° "' ] 


[mas yr^^] 


[mas yr^ 


14380 


1 33 21.800 


+60 40 12.25 


-0.29 


+0.74 


14561 


1 33 33.908 


+60 42 59.62 


-1.00 


-1.74 


14503 


1 33 30.184 


+60 39 19.62 


-0.68 


-0.53 


14411 


1 33 24.347 


+60 39 44.88 


+0.10 


+2.95 


14288 


1 33 16.193 


+60 39 43.50 


-0.51 


+1.27 


13899 


1 32 46.876 


+60 31 02.94 


+9.07 


+6.72 


13894 


1 32 46.417 


+60 46 26.62 


+0.04 


-1.44 


14223 


1 33 10.977 


+60 39 31.50 


+0.16 


-0.18 


14341 


1 33 19.427 


+60 40 47.50 


-0.25 


-1.30 


14928 


1 34 02.261 


+60 40 00.25 


-0.87 


-0.55 



According to the method presented by Sanders ( |1971[ ), we 
fitted a sharp (for the members) and a wider spread (for the 
field stars) Gaussian distribution to the distribution of the stars 
in the vector point plot diagram. We computed the parameters 
of the two distributions with a maximum likelihood method. 
From the values of the distribution at the location of the stars 
in the diagram we derived the membership probabilities. Due 
to the small difference between the centres of the two distribu- 
tions the stars with proper motions far away from the maximum 
are clearly identified as non-members, whereas it is difficult to 
decide which objects of the central region do belong to the clus- 
ter and which not. This is represented in the histogram for the 
membership probabilities (Fig.^ with a clear peak at P ~ 
and a much less pronounced increase towards P = I. Fig. |] 
shows the positions of the member and non-member stars. 

With this method 77 stars were classified to be members of 
NGC 581 with a probabiHty of at least 0.8, while 151 objects 
show a lower membership probabihty. The average proper mo- 
tions for cluster stars are: 

fiaCosS = (0.11 ± 1.07) mas yr^"'^ 
Us = (-0.95 ±1.24) mas yr^i 

and for the field stars: 

fiaCos6 = (1.14 ± 4.46) mas yr"^ 




/igCos d[mas/yr] 

Fig. 3. Vector point plot diagram of the stars in the region of 
NGC 581. Note that it is difficult to distinguish between the 
cluster proper motions and the centre of the field stars' proper 
motions. The stars with a membership probability of less than 
0.8 are indicated by small, the others by larger dots. The width 
of the distribution of the stars with a high membership prob- 
ability is of the order of 1 mas yr^^ which coincides with the 
standard deviation of the proper motion of a single star 




Fig. 4. Histogram of the membership probability for the stars 
from Fig. || 
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60.8 



60.6 - 




a[h] 



Fig. 5. Diagram of the {a,d) positions of the cluster mem- 
ber and non-member stars due to the proper motions. The non- 
member stars are indicated by small, member stars by larger 
dots. Note that the field stars are distributed homogeneously 
over the whole area, however, the member stars are located in 
a larger than expected area. This point is discussed in Sect. K2 



'0.91 ±3.77) mas yr^ 



Although the values are very similar, the standard deviations 
show a large difference between field and cluster stars. We 
should remark that the "true" members might be more concen- 
trated towards the cluster centre, but because of the problems 
in dividing between cluster and field stars, some non-members 
might be taken for member stars and vice versa. 

3.3. Tmmpler 1 

Near the edge of the photographic plates we found the open 
star cluster Trumpler 1 for which we attempted a proper motion 
study, too, as a by-product of our work. 

In a region of 0.15° around the centre of the cluster at 
02000 = 01*^35^.7, (^2000 = +61° 17', we detected only 64 
stars on the base of the photographic plates. CCD photometry 
is not available. A membership analysis with the same methods 
as described above resulted in average proper motions of 

/iaC0S(5 = (0.59 ± 1.71) mas yr^^ 
^ls = (-2.68 ±1.19) mas yr^i 

for the cluster members and 

/iQ,cos(5 ~ (1.41 ± 5.73) mas yr^^ 

= (-4.59 ±3.93) mas yr"^ 



/igCos d[mas/yr] 

Fig. 6. Vector point plot diagram of the stars in the region of 
the open star cluster Trumpler 1 . Again, the stars with a mem- 
bership probability of less than 0.8 are indicated by small, the 
others by large dots. Be aware that this proper motion study 
is of less reliability than the previous one of NGC 581, as it is 
based only on photographic plates, and the cluster is located 
close to the edge of the plates 



for the field stars. We present a vector point plot diagram in 
Fig. I 

It should be mentioned that these proper motions are less 
accurate than our results for NGC 58 1 , as Trumpler 1 is located 
close to the edge of the plates and we did not include CCD po- 
sitions in our analysis. However, the results indicate that Trum- 
pler 1 and NGC 581 have not only r oughly the same ages and 
distances (see Phelps & Janes 1993| ), but — within the errors 
— the same absolute proper motions. Therefore, the two ob- 
jects may be taken as a candidate for a binary star cluster (see 
e.g. Subramaniam et al. 1995 ) as they are known for the Mag- 
ellanic Clouds (Dieball & Grebel |1998|). 



4. Analysis of the colour magnitude diagram 

4.1. Field star subtraction 

The CMD derived in Sect. ^ is contaminated with field stars. 
Before trying to analyse the CMD it is essential to subtract 
these stars to emphasise the features which belong to NGC 581 
and to prepare the data for the determination of the cluster IMF. 

In the magnitude range covered by the proper motion study 
it is possible to distinguish between field and cluster stars using 
the information of the membership probabilities from Sect. ^ 
The proper motion study is complete down to = 14 mag so 
that until this point, a membership probability of 0.8 or higher 
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Fig. 7. Luminosity functions down to y = 14 mag according 
to field star subtraction with proper motions (solid line) and sta- 
tistical subtraction (dotted line). Statistical field star subtraction 
for the brighter stars would have lead to 9 additional cluster 
members and therefore to a shallower IMF 



is considered a suitable criterion for the definition of the mem- 
bers of NGC 581. 

Below this limit, the star numbers are high enough to jus- 
tify a statistical field star subtraction. For this, the CCD field 
was divided into two regions of the same area of 1.3 • 10^ pix- 
els: a circular region with the star cluster in its centre (radius: 
654 pix. or 18') and a ring encircling this region (outer radius: 
925 pix. or 26')- Under these assumptions, all cluster members 
are surely located in the central area. For both parts separate 
CMDs were computed. Each of them was divided into rectan- 
gular bins with a length of 0.5 mag in magnitude and 0.1 mag 
in colour. (Variation of the bin sizes did not affect the results.) 
The numbers of stars in the field CMD cells were determined 
and as many stars of the corresponding cell of the inner region 
CMD were randomly chosen and removed. Assuming a homo- 
geneous distribution of field stars over the whole field of view, 
the resulting CMD, which is presented in Fig. H, represents only 
the cluster members. 

To illustrate the advantage of our method. Fig. ^ shows 
the luminosity functions obtained with the proper motions and 
the statistical field star subtraction, respectively, in the range 
of the photographic plates (V < 14 mag). It can be seen that 
the proper motion analysis leads to less member stars than the 
statistical subtraction. 



4.2. Morphology of the cluster CMD and age determination 

The brightest stars in the field of NGC 581 (marked with let- 
ters A to C in Fig. [l|) were saturated even in the shortest expo- 
sures so that no direct photometric information was available 
for them. However, the three stars were the subject of previous 
photoelectric studies (Hoag et al. 1961 , Purgathofer 1964 ) so 
that their V magnitudes and B — V colours could be added to 
our CMD manually. Star A was too bright to be measured on 
the photographic plates with a sufficient accuracy so that we 
could not derive its proper motion, either However, since star 
A is a HIPPARCOS stai- (ESA |1997|, stai- no. 7232), we could 



get its proper motion from this source. Calculating its mem- 
bership probability led to a value of 0.89. Star B has a proper 
motion far away from the one of the cluster so that it can be 
considered a non-member, while star C belongs to NGC 581 
with a probability of 0.90. 

After field star subtraction, we found three red stars with 
B ~V k2 mag from V « 10.8 mag to y « 12.4 mag. Com- 
pared with the well populated main sequence in the same mag- 
nitude range, it seems obvious that these objects cannot be clus- 
ter members. We assume these objects to be field stars which 
are classified as member stars as a consequence of the small 
difference between field and cluster proper motions. 

The secondary main sequence and giant branch like struc- 
tures completely vanished within the range of the proper mo- 
tions. Below V — lA mag, we still find some stars remaining 
in the region of these features. Since the vector point plot di- 
agram does not give any evidence for a second cluster in the 
same line of sight, we assume that these stars belong to the 
field star population(s). It is very unlikely, too, that the remain- 
ing stars are pre main sequence cluster members, because at the 
age of NGC 581 (see the following paragraph) this kind of ob- 
jects should be much closer to the main sequence than the stars 
in the CMD are (Iben 1965). Therefore we assume these stars 
to be remnants of the field star subtraction due to the imperfect 
statistics of the sample. 

We fitted isochrones of the Geneva group (Schaller et al. 



1992 ) to the resulting CMD from which we derived the distance 



modulus, reddening, age, and metalicity of NGC 581. The best 
fitting isochrone is plotted into the cleaned CMD of Fig. ^ 
The parameters of the selected isochrone are sho wn in Table 
H Isochrones of the Padua group (Bertelli et al. 1994 ) were 
fitted to the CMD for comparison and led to the same set of 
parameters. 

4.3. Completeness correction 

To obtain comprehensive luminosity and initial mass functions, 
the data have to be corrected for completeness. As crowding is 
not a problem in our images, the completeness in the field and 
cluster regions are the same so that it is not necessary to correct 
for completeness before the field star subtraction. 

With artificial star experiments using the DAOPHOT II rou- 
tine addstar we derived the completeness function shown in 
Fig. ^. As the V images reach down to much fainter magni- 
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Fig. 8. Colour magnitude diagram of all members of NGC 581 
as determined with the proper motions (V < 14 mag) and the 
statistical field star subtraction (V^ > 14 mag). The parameters 
of the isochrone plotted in the diagram are listed in Table ||. 
Note that the brightest star of Fig. || is not a member of the 
cluster 

Table 8. Parameters of NGC 581 as derived from isochrone 
fitting to the colour magnitude diagram. The errors of distance 
modulus and reddening are estimated, the uncertainty of the 
age is derived from the comparison with the isochrones neigh- 
bouring the selected one 



distance modulus 


(m - M)o 


= 12.3 ±0.1 mag 


i.e. distance 


r 


= 2884 ± 130 pc 


reddening 


Eb-v 


= 0.41 ± 0.02 mag 


age 


\ogt 


= 7.2 ±0.1 


i.e. 


t 


= 16 ± 4 Myr 


metalicity 


Z 


= 0.02 



tudes, these experiments were only performed with the B ex- 
posures. Down X.0 B — 17.25 mag the sample is more than 
80% complete with a sharp drop afterwards to almost 0% at 
B = 19 mag. Figure || shows that the completeness level of 
60% is reached around B = 1% mag, so that we assume this 
value to be a reasonable cut-off for our studies. With a main 
sequence star colour of B — V = 0.8 mag for B = 18 mag, 
this corresponds to a limiting magnitude ofV = 17.2 mag. 

4.4. Initial mass function 

With the stars of the cleaned CMD we determined the lumi- 
nosity (Fig. |l^) and initial mass functions of NGC 581 after 
deleting all objects outside a A{B — V) — 0.3 mag wide 



Fig. 9. Diagram of the completeness of our data as a function 
of B magnitudes. Down to B = 17.25 mag, the sample is more 
than 90% complete, down to one magnitude fainter, we detect 
a completeness of more than 50%. After this point, the com- 
pleteness sharply drops until no more stars are being detected 
at B ~ 19.5 mag. For our further investigations we skip all 
stars fainter than i? = 18 mag which corresponds to a com- 
pleteness level of higher than 60% 

strip around the main sequence and above the turnover of the 
isochrone around V = 10 mag from the CMD and applying the 
completeness correction described in Sect. 4.3 As mentioned 
in Sect. 3.2, some field stars will have stayed in the sample 
while a few cluster stars might have been rejected. However, we 
assume these two effects do not influence the IMF of NGC 581. 

From the initial stellar masses given in the Geneva 
isochrone data, we computed a mass-luminosity relation rep- 



resented by a 6*^^ order polynomial: 

6 

m[MQ]^Y,d,-V'[mag] 



with the following parameters: 



(7) 



do 
di 
d2 
d3 
di 
db 
de 



-523.815 
±221.123 
-36.347 
±3.055 
-0.140 
±0.003 
-0.0000327. 

A 5'''^ order polynomial still did not fit the faint part of the 
mass-luminosity relation well enough. With Eq. (Q) we derived 
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Fig. 10. Luminosity function of NGC 581. The histogram ob- 
tained after completeness correction is printed with solid, the 
uncorrected one with dotted lines 



the initial masses of the stars classified as objects of NGC 581 
on the base of their V magnitudes. V was preferred compared 
with the B values as their photometric errors are smaller at 
equal magnitudes. 

We included all stars with a mass of higher than m = 
\A\Mq (corresponding to the previously mentioned value of 
V = 17.2 mag or log to — 0.15). In this range, the complete- 
ness is 60% or higher. 198 stars were selected by this criterion. 
The stars were divided in A log to = 0.1 wide bins. The star 
numbers were corrected using the results of the completeness 
calculation, resulting in 247 "virtual" stars from 9.45Mq down 
to 1.41Mo. The resulting histogram is presented in Fig. [Til . 
The single star with a mass around 9.45Af0 (log to = 0.98) 
was not taken into account for the IMF determination, since 
this single object at one end of the histogram might heavily in- 
fluence the IMF slope. The slope of the IMF was determined 
by linear regression to the histogram. 

We obtained T = -1.80 ±0.19 from the stars with a 
mass from IAIMq to 8.7OM0 (An experiment including the 
9.45M0 star led to F = -2.11 ± 0.23 and therefore to a dif- 
ference in the exponent of more than 0.3!). The resulting IMF 
is shown in Fig. |ll[ 

5. Summary and discussion 

NGC 581 is a young open star cluster with an age of 16 Myr. 
It is located at a distance of 2900 pc from the Sun. We derived 
proper motions of 228 stars in the region of the cluster down to 
14 mag. 77 of those can be considered members of the object. 



Fig. 11. Initial mass function of NGC 581. The solid line corre- 
sponds to the completeness corrected values, the dotted line to 
the original data. The IMF slope calculated with linear regres- 
sion to the histogram is determined to be F = —1.80 ± 0.19 

A study of the IMF of NGC 581 leads to a power law with a 
slope of r = -1.80. 

The results of the photometry mainly coincide with the 
findings of Phelps & Janes (1993 1, however, we do not find ev- 
idence for a star formation over a period of as long as 10 Myr: 
While Phelps & Janes claim the necessity of two isochrones of 
different ages to fit both the blue and red bright stars in their 
CMD, the Geneva \ogt = 7.2 isochrone fits all features of our 
CMD sufficiendy well. 

The IMF slope might be slightly higher t han Sc alo's ( |1998| ) 
figure for this mass range, as Phelps & Janes ( 1 994 ) find a steep 
slope for NGC 581, too. Their value of F = -1.78 marks 
the steepest IMF of their entire study of open star clusters. 
Although their IMF is based on a photometry with a deeper 
limiting magnitude, our determination has its advantages, too: 
Phelps & Janes only used statistical field star subtraction and no 
information of proper motions, which may have great impact 
on the high mass range and therefore might affect the slope of 
the IMF. Furthermore, they only took one single region in the 
vicinity of their star clusters for the field star subtraction, while 
we were able to use the surrounding of the cluster itself. And 
finally, their field of view is so small that they probably did not 
cover the entire cluster which, assuming the presence of mass 
segregation, might have flattened their IMF due to a lack of low 
mass stars from the outer regions of the cluster 

If we would not have based part of the membership deter- 
mination on the proper motions, we would have ended up with 
9 more stars in the magnitude range of < 14 mag (see Fig. 
0). Since the total star number in this region is low, this would 



have decreased the IMF the slope to T = 1.73 ± 0.28. Within 
the errors, this result would still match with our value above, 
however, the shape of this IMF is less well represented by a 
power law, which is expressed by the higher error Therefore 
we claim that adding the information of the field star subtrac- 
tion did improve the reliabiUty of our IMF study. 
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